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■ Abstract 

o : 

O ; 

\ Recently Garcia Bellido et. al. have proposed that electroweak baryogenesis 

may occur at the end of inflation, in a scenario where the reheat temperature 
' is too low for electroweak symmetry restoration. 1 show why the scenario is 

I ■ difficult to test reliably by classical field techniques on the lattice. 

(N : 

^ ; I. INTRODUCTION 

> . 

O ■ Inflation is a plausible and popular explanation for how the initial cosmological epoch 

■ produces a universe with such startling size, flatness, and homogeneity ||1|]. However, inflation 
O ■ makes even more puzzling the other remarkable feature of the universe-that it contains a 
^ \ macroscopic but relatively small net abundance of baryonic matter (approximately 5 baryons 
^ ; per 10^° photons [|). 

p • Recently Garcia-Bellido, Grigoriev, Kusenko, and Shaposhnikov have proposed a way 

p • that certain inflationary scenarios may be able to explain baryogenesis (the origin of the 
D . baryon abundance) as well [Q. During inflation, most of the energy density in the universe 
\ is in the potential energy of a scalar fleld, the inflaton. It has recently been understood that 
\ inflation can end, and the energy density stored in the inflaton fleld can be converted into 
^ ■ a thermal bath, much more abruptly than had previously been thought possible, a process 
' called "preheating" [Q. In the scenario of Garcia-Bellido et.al., the baryons are created 
during this far from equilibrium process. 

Baryon number is not conserved in the standard model f^. The violation arises from 
nonperturbative physics of the SU(2) weak gauge flelds (W and Z bosons). When occupation 
numbers of infrared flelds become large, nonperturbative physics can become efficient. This 
happens at high temperature and can also happen in other high excitation situations. In 
equilibrium, it occurs when there is no Higgs fleld condensate,[] which requires a temperature 
T > lOOGeV. At lower temperatures there is a Higgs fleld condensate and baryon number 
violation is exponentially slow 0. In the Garcia-Bellido scenario, the energy density in the 
inflaton passes flrst into very infrared fleld modes, and baryon number is readily violated; 



^The notion of a Higgs field condensate is a perturbative one, and should be used with some caution. 
Nonperturbatively speaking there is no qualitative distinction between symmetry broken and restored phases 
1^ and there can either be a phase transition or an analytic crossover between them, depending on coupling 
parameters 0. 



but when the fields fully thermalize the temperature is low enough that there is a Higgs field 
condensate and no further baryon number violation occurs. 

The Garcia-Bellido scenario involves nonlinear, nonperturbative physics, and its quan- 
titative study is difficult to conduct analytically. For preheating in general, classical field 
techniques have proven useful P,p!0|. Baryon number violation in equilibrium has also been 
treated accurately by classical field techniques [p!lH T5[|, and recently these techniques have 
been applied to the Garcia-Bellido scenario, both in a 1+1 dimensional toy model P,p!B| and 
in a more realistic 3+1 dimensional setting |[17| . 

It sounds natural to expect classical field techniques to work well for the study of the 
Garcia-Bellido scenario. This note will argue that there are serious complications, because 
in the context of 3+1 dimensional Yang-Mills theory at realistic coupling, the techniques as 
they exist to date contain spurious physics which can lead to "fake" early thermalization 
and baryon number violation. At best, classical field techniques will probably have to be 
modified and used with care in this context; at worst they may not be useful at all. 



II. CLASSICAL YANG-MILLS AND THE LATTICE 

A. Classical field approach 

First I summarize how classical field techniques are applied to preheating. The early 
stage of preheating can be understood analytically. When the infiaton field condensate 
begins to oscillate about its potential minimum, certain field modes coupled to it have 
their field amplitudes grow exponentially due to parametric resonance @]. In l|,|l3 it is 



shown that classical fields, initially populated with "quantum vacuum initial conditions," 
give quantitatively the same behavior. Here "quantum vacuum initial conditions" means 
Gaussian random initial conditions in k space, with mean squared excitation equal to the 
vacuum zero-point excitations of the relevant modes, {h = 1) 
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with Q and P real canonical field and momentum variables. When parametric amplification 
has progressed the physics becomes nonlinear; but the occupation numbers of the relevant 
field modes are large by then, and the classical field approximation is valid. The classical 
field theory can be discretized on a lattice and evolved in real time by standard algorithms. 
This takes fully into account the nonlinear interactions of different field modes. Provided 
all of the interesting physics involves modes well to the infrared of the lattice spacing scale, 
the discretization should not disturb the physics of interest. 

Such "quantum" initial conditions will not be preserved by classical field evolution; noth- 
ing prevents the energy associated with the "zero point" fiuctuations, which have been 
turned into excitations of classical fields, from moving between Fourier components of the 
field. They will not do so, or will do so very slowly, if the coupling constant is small. But 
at finite coupling, even without an oscillating scalar (infiaton) condensate, this is not so. 
Instead, the lattice, classical field evolution will eventually approach the classical thermal 
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distribution for the lattice system in question, which at weak couphng has0 

2 ^2 

{Q ( ^)) thermal latt — T, (-P ( ^ )) thermal latt = T , (2.2) 

/C^ _|_ YYi-i 

for some T ~ 1/a. (By fc^ I mean the lattice dispersion relation, a^k"^ = 4 sin^(/cja/2).) 
Note that the distribution is not Lorentz invariant-Lorentz invariance is broken by the lattice 
discretization-and has sma//er field excitations than the vacuum ones at the largest k"^, but 
much larger field excitations in the infrared. 

B. Yang-Mills theory: Coulomb gauge correlators 

This is potentially problematic for the case at hand because the very infrared physics is 
the physics of interest. Eventually, the energy in the UV modes will cascade into the IR 
modes in a manner impossible in the actual quantum theory of interest. The question is: 
On what time scale does this occur? We should test whether this takes place on a time scale 
short enough to disturb the physics which is of interest to the simulations. The theory of 
interest is SU(2) Higgs theory in 3+1 dimensions, at realistic coupling. A simple way to test 
the viability of the classical technique is to examine the case where there is no condensate, to 
see how the "vacuum" fields evolve classically in isolation. If the time scale for poor behavior 
is shorter than the time scale for "interesting" behavior when the djTiamics of interest is 
present, there is a problem. 

We need some measurables for our study, which can indicate whether the cascade of 
energy into the IR has occurred; I choose to study the equal time two-point function in 
Coulomb gauge, as a function of Coulomb gauge is not expected to give sensible results 
for unequal time correlators because the gauge fixing procedure treats each time slice sep- 
arately, except for a single global time dependent gauge freedom. However, for equal time 
correlators it should give sensible results, although it is unclear exactly how to interpret the 
most infrared modes. For wave numbers where perturbation theory is useful, we can inter- 
pret the values of the Coulomb gauge fixed correlation functions in terms of a A; dependent 
"effective temperature," similar to what is done in scalar field studies of inflation. We see 
from Eq. ( p.2[ ) that the effective k dependent temperature is 

Teff(fc,A) = ^Ck' + m'){AUk)) , Teff(A:,E) = ^{Efik)) , (2.3) 

where the factor of 1/2 accounts for the two transverse modes summed over in the index 
trace; the longitudinal mode of A is fixed to zero by the Coulomb gauge fixing procedure. 



^Thcrc arc corrections to (Q^) of relative size AT/Vfe^ + m^, due to interactions. However, if P enters 
the Hamiltonian only quadratically and is not involved in constraints (such as Gauss' Law for E in gauge 
theories) then the expression shown is exact and defines T. 

We fix Coulomb gauge with the definition of Mandula and Ogilvie ]l8[: b ut there are algorithmic differences 
appropriate for the real-time setting. Our algorithm is discussed iiT~|l3|. The lattice implementation and 
classical field evolution algorithm follow Ambj0rn et. al. [ pT| . 
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FIG. 1. Equal time correlators in Fourier space, for classical lattice SU(2) theory at realistic coupling 
(7w = 0.64, evaluated in Coulomb gauge, at a series of times after "vacuum" initial conditions are drawn as 
described in the text. Top panels show (fc^v4^(fc, t)), bottom panels show {E'^{k, t)); up to interaction effects, 
the results can be interpreted as "fc dependent temperatures." The right panels show a close up of the small 

region, and earlier times. In equilibrium at zero coupling each curve would be flat; the equilibrium E 
field curve is not flat because of Gauss' Law, and the A field curve is not flat because of interactions. 
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and the longitudinal mode of E is set zero, to leading order in perturbation theory, by the 
Gauss constraint. The above procedure also defines T^s even where perturbation theory 
does not apply, but it is not clear how to interpret the result. Due to interactions, even in 
equilibrium Tefj(A;) will not be constant for all k, nor will it be equal for A and E fields; but 
wherever the behavior is close to free field behavior, Tes{k) will be approximately fiat and 
approximately the same for A and E fields. 

Since the cascade of energy, in the gauge sector, should be at least as fast in the full 
theory as in just Yang-Mills theory, and since Yang-Mills theory has fewer parameters, I 
study just Yang-Mills theory. Hence, = 0. The initial conditions are chosen as follows. 
All field and momentum variables are chosen as in Eq. ( |2.1| ), then the E fields are projected 
to the Gauss constraint surface; the presence of longitudinal A field excitations in this choice 
of initial conditions is irrelevant because fixing to Coulomb gauge will remove them. This 



is the same procedure for choosing initial conditions as was used by Rajantie et. al. [|17 
Note that there is only one length scale, a, in the problem, and the gauge coupling is a 
independent, so all behavior scales as a is changed; the behavior of a mode of wave number 
k depends only on ak. 

Figure shows how the field excitations evolve with time when we choose initial condi- 
tions as in Eq. ( |2.1D at gauge coupling g^i; = 0.64, on a 64^ lattice. The striking feature in 
the figure is that the infrared excitations (fc^a^ < 1/3) are populated nearly to their thermal 
value very early after the evolution begins; they are largely populated by t = 8a and fully 
populated at t = 16a. Then, over a much longer time scale, the ultraviolet modes thermalize 
among themselves. 



C. Diffusion of Chern-Simons number 



A possible objection to this gauge theory study is that it focused on gauge fixed ob- 
servables. One should always treat gauge fixed measurements with care, so I also study 
baryon number violation directly. The chiral anomaly relates the baryon number to the 
Chern-Simons number, 



-NniA,) = NcsiA, 



9' 



A'=A, 



dr d X 



Bt{A')D,A 



(2.4) 



where r parameterizes a path through the space of 3-dimensional gauge fields from the 
vacuum to the gauge field of interest. As defined here A^cs modulo 1 is gauge invariant; 
the integer part depends on the path choice. Choosing the path to coincide with the real 
time field evolution, the ensemble average of A'^S' (^cs)) ^^^^ grow linearly in time when 
baryon number violation is active. Since A"^g also has a large "noise" contribution from UV 



excitations [0 it is more convenient to measure A^cs not of the actual configurations in a 



field evolution but of "cooled" copies; this is discussed in |]T4| , and extensively in [jl9| , which 
also describes the algorithm used here. 

Figure shows (A^qs) as a function of time, averaged over about 1000 trajectories with 
independent initial conditions, on a 32^ lattice and using a cooling depth (defined in |]T9| ) 
of 7.5a^. It also shows the average of the square of the integer part of A^cs (the stars, which 
are based on a subset of about 400 of the trajectories, and therefore have correlated and 
somewhat larger errors). The integer part is defined as the difference between evaluating 
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FIG. 2. Diffusion of Chern-Simons number, starting from "vacuum quantum" initial conditions, in a 
32'^ box at = 0.64. The band of error bars show the mean square value of A^cs as defined in the text; 
the stars indicate the mean squared values of the topological, integer part of iVcs. Efficient baryon number 
violation begins at i ~ 10a. Data at all t values are from the same set of time evolutions, so the errors are 
highly correlated; the integer measurement was performed on a subset of the data used for the band of error 
bars, and have errors about 50% larger. 



Eq. ( p.4|) for A^cs(^) using the time history as the path, and using the cooling path from 
Ai{t) to the vacuum as the path. It is a topological number which directly measures whether 
genuine topology change has occurred. If the vacuum state were preserved by the time 
evolution, we would expect (iV^g) to fluctuate about an average close to zero, and the 
topological number would remain strictly zero. Instead, the flgure shows that there is a 
brief initial period where this occurs; then, by t = 10a, Nqs begins to diffuse at roughly the 
same rate it does after full thermalization. This means that, at the physical value of the 
coupling = 0.64, the IR gauge flelds relevant for baryon number violation thermalize by 
about if: = 10a, and baryon number proceeds from there. This effect is entirely an artifact 
of using classical lattice evolution with "quantum zero point" initial conditions rather than 
treating the actual quantum theory, where baryon number, in vacuum and at this coupling, 
is exponentially suppressed with a large exponent. The time scale for this "artifact" baryon 
number violation to begin quite short; in particular it is shorter than the time scale on 
which Rajantie et.al. report nontrivial baryon number violating physics in the Garcia-Bellido 
]r^], by almost a factor of 3. This brings seriously into question the reliability of 



scenario 
their results. 
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FIG. 3. Time development of equal time (j) correlator in Fourier space for two component scalar field 
theory. Thermalization is slow and does not proceed first for the infrared modes, as it does in the gauge 
theory case. 

D. Comparison to the scalar theory case 

For comparison, and because it bears on the classical field theory techniques applied in 
preheating P,p!0 



I have also analyzed a two component scalar field theory with Lagrangian 



2 \ 



(2.5) 



This theory behaves very differently from the gauge theory; indeed, even at quite large 
coupling the time scale for thermalization turns out to be very long. In Figure |^ I show how 
the power spectrum (which is now a gauge invariant quantity) evolves at the quite large 
coupling value of A = 1, and bare mass squared a^m^ = — .82A, chosen to approximately 
cancel a "tadpole" contribution so the initial behavior is close to that of a massless theory. 
The figure shows Tes{k, P), the effective temperature derived from the momentum degrees of 
freedom, only; in thermal equilibrium this will be a flat line (up to statistical fluctuations) at 
any value of the coupling. The flgure shows that, in the scalar field case, the thermalization 
is very slow, even at A = 1, and it proceeds across all A;^ values at about the same rate. 
Since the preheating literature generally deals with much weaker couplings, this suggests 
that classical lattice techniques can be reliable in the preheating context, unless the time 
scale of the interesting physics is extremely long. 

III. DISCUSSION 



We have shown that classical Yang-Mills theory, regulated on a lattice and given "quan- 
tum vacuum" initial conditions, shows very rapid heating of the infrared modes, and then 
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more gradual thermalization between all modes. The time scale for heating of the infrared, 
for SU(2) at realistic coupling = 0.64, is of order 10 lattice spacings of time (independent 
of the lattice spacing). This is probably faster than any purely infrared, nonequilibrium dy- 
namics of interest will develop, which makes it difficult to study baryogenesis at preheating 
without the results being contaminated by artifacts of the lattice technique. 

The physical reason for the rapid transfer of energy between hard and soft degrees of 
freedom can best be understood in the language of plasma physics. The excitations es- 
tablished for ultraviolet degrees of freedom, intended to simulate their quantum zero point 
fluctuations, propagate nearly ballistically. However, because the theory is nonabelian, they 
carry nonabelian "charge," and constitute a plasma in which the IR fields evolve. The 
plasma degrees of freedom (UV lattice modes) move in the background of the IR modes of 
interest, and infiuence their evolution rather as a plasma modifies the evolution of infrared 
electromagnetic fields. It has long been appreciated that such a plasma strongly modifies 
the dispersion relations of the soft excitations (Debye screening, plasma oscillations, etc.), 
and efficiently changes the amplitude of soft excitations via Landau damping. In the non- 
abelian context, for quantum, equilibrium situations, this physics is contained in the "hard 
thermal loops" of Braaten and Pisarski ||2^. An extension to the classical, lattice setting 
has been addressed by Arnold The damping away of large gauge field excitations is 
efficient below the scale k ^ g/a {in a. thermal plasma the scale would be ~ gT). However, 
damping is a two way street; the time scale for a large IR field to be damped away is also 
the time scale for UV couplings to populate that IR field, at a "temperature" corresponding 
to some average over the effective temperatures of the UV degrees of freedom responsible 
for the damping. In the simulations of this paper, no large initial field exists; so we only see 
the populating of the IR modes, to an "effective temperature" of order an average over the 
temperatures implied by the excitation levels of the UV modes. This is approximately what 
we see in Figure |I|; the IR modes with k < gT are quickly excited to an effective temperature 
comparable to the effective temperature of the UV modes; then over a longer time scale, 
less infrared modes also thermalize. 

Are these lattice artifact physics effects relevant for the study of the Garcia-Bellido sce- 
nario, where besides the "vacuum" excitation I have addressed, there is also large, coherent 
initial infrared field physics? I will argue that they are, and that they seriously complicate 
the analysis of the Garcia-Bellido proposal by lattice means. 

The interesting physics in the Garcia-Bellido scenario is infrared, large occupation num- 
ber physics. We must make the length scales associated with that physics (say, the wave- 
length of modes driven on resonance) much larger than the lattice spacing, or else the 
interesting physics will directly be contaminated by lattice effects. Therefore the time scale 
on which the nontrivial IR dynamics is expected to occur is quite generically long compared 
to the lattice spacing a. For instance, the interesting dynamics in the simulations of Rajantie 



et. al. |]T^ occurs at t ~ 30a after the beginning of the simulation. Therefore we can expect 



that the excitation of the IR gauge field modes by their coupling to the UV should already 
have taken place before the interesting IR physics is complete. 

If the actual (quantum, nonequilibrium) dynamics of interest fails to excite large IR 
gauge fields, then in the simulations we will nevertheless see large IR gauge fields excited, 
by the coupling to the UV described here. There is therefore clearly the possibility of a "false 
positive" indication of baryon number violation (if for instance Higgs symmetry is briefiy 
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restored, but in the actual dynamics the gauge fields would be too cold to violate baryon 
number). Perhaps less obviously, it is also possible that the coupling to the UV lattice 
modes will create "false negative" results where baryon number violation actually should 
occur. This could happen if the true (quantum, nonequilibrium) dynamics actually excites 
the IR gauge fields very efficiently, to an effective temperature higher than what the UV 
modes would supply. In this case, the UV modes represent an efficient absorber of the IR 
gauge field excitation energy, via Landau damping. In the real (vacuum, quantum) theory, 
energy loss to the UV should not be too efficient, and should occur mainly by a cascade. On 
the lattice, the UV lattice modes can directly Landau damp away large IR gauge fields on a 
time scale ~ 10a, short compared to the (expected) dynamic timescale for the nontrivial IR 
dynamics. Note that the UV modes also change the dispersion of the lattice modes, which 
means that the evolution may be wrong even if energy is not transferred to the UV. One very 
basic way of checking for some of these problems is to test for lattice spacing dependence in 
the results. This was not done in |]T7| . 

It is possible that there are ways of evading the problems discussed here. For instance, 
one could choose initial conditions in which only excitations with /c^ below some cutoff 
were excited. But it is not obvious that this treatment will reproduce the correct (quantum) 
treatment either. The physics of interest probably involves energy moving from a condensate 
into IR modes, and then cascading into the UV. Will this cascade be incorrectly described if 
the UV modes start out with no excitation? Does failing to quantize the UV modes already 
ensure that the description cannot be correct? The burden of proof clearly lies with the 
practitioner. It would be necessary but not sufficient to demonstrate that all simulation 
results show weak dependence on the choice of fc^ cutoff. 

The behavior we have seen should not he expected in 1+1 dimensional, abelian studies 
of baryogenesis at preheating, such as those of Garcia-Bellido et. al. P,p!^. This is a matter 



of dimensionality. In 3+1 dimensions, thermal self-energy corrections are UV divergent in 
a classical theory, with the divergence cut off in nature at ~ T by quantum mechanics. 
However, in 1+1 dimensions thermal self-energies are UV finite. Therefore the studies of 
the 1+1 dimensional abelian analog theory in |^,16| are deceptive. The behavior I describe 
also appears to play little role in abelian Higgs theory, at least at the unrealistically small 
coupling considered in p2|. I should also emphasize that my results do not mean that 
previous studies of infiationary preheating are incorrect. As Figure ^ shows, the time scale 
for UV energy to cascade to the IR in a scalar theory, even at the very large coupling of 
A = 1, is very long, thousands of lattice units of time. The cascade time is expected to 
grow at weak coupling as A~^. The largest value of A used in |1^ was about two orders of 
magnitude smaller than that in Fig. |^, so the time scale for the cascade to occur in their 
work would be at least lO'' lattice units, much longer than any time scale they considered. 
Physically this difference is related to the fact that the hard thermal loop for a gauge field 
contains an imaginary (Landau damping) part, while that for a scalar field does not. 
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